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Prefacio 


En los últimos años, el uso de MATLAB9! como lenguaje de programación y en- 
torno de computación científica se ha generalizado tanto en el ámbito académico como 
en la industria. Su flexibilidad en la implementación de código y las herramientas de 
visualización de que dispone le han convertido en un estándar para la resolución de 
problemas no triviales provenientes de diversas ramas de las ciencias aplicadas. 

Desde el punto de vista docente, MATLAB constituye una herramienta muy versátil 
para iniciar a los alumnos de cualquier disciplina científica en el campo de la simu- 
lación numérica. En esencia, este texto tiene como finalidad establecer los fundamentos 
generales de MATLAB para que, posteriormente, cada alumno pueda profundizar en 
las aplicaciones específicas que más le interesen. 

Los temas tratados a lo largo del texto engloban aspectos básicos de funcionamiento 
del programa, que van desde la definición de vectores y matrices hasta la implementación 
de funciones por parte del usuario a través de los denominados M-files. Otras cuestiones 
más específicas incluyen la descripción del módulo de cálculo simbólico (Symbolic Math 
Toolbox) y una exposición detallada de las herramientas para realizar representaciones 
eráficas en dos y tres dimensiones. El texto incorpora ejemplos y ejercicios de apli- 
cación de los comandos introducidos. No se asumen conocimientos previos de MATLAB, 
aunque es recomendable que el lector conozca los principios básicos de la programación 
para un mejor aprovechamiento de los contenidos descritos. 

Finalmente, se ha incluido una sección de referencias bibliográficas que contiene, 
en Opinión de los autores, algunos de los textos más didácticos para el aprendizaje de 
MATLAB. De todos ellos, cabe destacar, por su claridad expositiva y el espectro de 
contenidos que abarcan, los libros Numerical computing with MATLAB, de Cleve B. 
Moler, creador de la primera versión del programa, y MATLAB quide, de Desmond J. 
Higham y Nicholas J. Higham. Algunas referencias de libre divulgación completan la 
bibliografía básica relativa a MATLAB y sus diversas aplicaciones. 


Madrid, febrero de 2014 Andrés Arrarás 
Danilo Magistrali 
Laura Portero 


IMATLAB es una marca registrada de The MathWorks, Inc. 


l. ¿Qué es MATLAB? 


MATLAB, acrónimo de Matrix Laboratory, es un sistema interactivo para el cálculo 
científico y un lenguaje de programación especialmente diseñado para la implementación 
de algoritmos numéricos. Su primera versión data de 1978 y fue escrita por el analista 
numérico Cleve B. Moler en lenguaje Fortran. Posteriormente traducida a C en 1984, 
su principal finalidad era proporcionar un acceso sencillo a las librerías LINPACK y 
EISPACK, que contenían algoritmos clásicos del álgebra lineal numérica. 

El perfeccionamiento de su interfaz de usuario y sus salidas gráficas permitió al 
programa adquirir una gran popularidad en el ámbito académico. Las aplicaciones de 
MATLAB se fueron extendiendo a diversas ramas de las ciencias aplicadas y la in- 
geniería, gracias al desarrollo de toolboxes: librerías escritas en el lenguaje propio de 
MATLAB y destinadas a resolver problemas provenientes de áreas específicas. 

En la actualidad, MATLAB se utiliza como herramienta de apoyo en campos tan 
diversos como la mecánica de fluidos computacional, la estadística, la simulación de 
sistemas dinámicos, el procesamiento de señales, la visión artificial, el diseño de sistemas 
de control o la biología computacional. 


2. Entorno de trabajo 


Al iniciar una sesión de trabajo en MATLAP?, el programa muestra por defecto las 
siguientes ventanas: 


o Command Window (ventana de comandos): es la encargada de recibir las instruc- 
ciones de MATLAB. El símbolo >> se denomina prompt e indica que el programa 
está preparado para aceptar órdenes. Para ejecutar un determinado comando, 
basta con pulsar la tecla Enter. Adicionalmente, las teclas $ y ] nos permiten 
recuperar entradas anteriores. 


o Command History (historial de comandos): su función es almacenar automática- 
mente las órdenes ejecutadas en la ventana de comandos. Mediante la copia y el 
pegado habitual o con el uso del ratón, es posible recuperar instrucciones tecleadas 
previamente. 


o Current Folder (carpeta o directorio actual): muestra los ficheros incluidos en 
la carpeta de trabajo del usuario. Cualquier fichero que se quiera ejecutar debe 
aparecer en la carpeta actual. 


o Workspace (espacio de trabajo): contiene la totalidad de variables definidas du- 
rante una sesión de trabajo. Haciendo doble click con el ratón sobre una variable, 
MATLAB nos la muestra en formato de hoja de cálculo, a través de un editor 
denominado Variables. 


2En la última versión disponible, MATLAB 8.2 (R2013b). 
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MATLAB trabaja por defecto en aritmética de coma flotante de doble precisión. 
No obstante, el usuario puede seleccionar el formato de salida mediante los comandos 
siguientes: 


O 


format short: es el formato por defecto y utiliza un máximo de tres dígitos para 
la parte entera y cuatro para la parte decimal; si el número que se quiere mostrar 
necesita más dígitos, hace uso de la notación científica (e.g., 0.4327e-002); 


format long: utiliza un máximo de dos dígitos para la parte entera y quince para 
la parte decimal; si el número que se quiere mostrar necesita más dígitos, hace 
uso de la notación científica (e.g., 3.141592653589793e+002); 


format short e: utiliza notación científica con cuatro dígitos decimales; 
format long e: utiliza notación científica con quince dígitos decimales; 


format compact: reduce el espacio en blanco existente entre las instrucciones 
introducidas por el usuario y los resultados proporcionados por MATLAB. 


Enumeramos a continuación algunos detalles adicionales que pueden resultar de 
interés: 


3. 


MATLAB distingue entre mayúsculas y minúsculas, de modo que, como veremos 
a continuación, las variables a y A son distintas; 


todos los comandos de MATLAB se escriben en minúsculas, con sus argumentos 
entre paréntesis y separados por comas; 


si añadimos punto y coma (;) después de una instrucción, ésta se ejecuta pero el 
resultado no se visualiza por pantalla; 


la orden c1c limpia la pantalla de comandos; 
la combinación de teclas Ctr1 + C hace que se interrumpa una ejecución; 


el comando diary nombre_fichero.txt permite grabar en un fichero de texto el 
trabajo realizado durante una sesión; con las instrucciones diary on y diary off 
activamos y desactivamos la grabación en el fichero, respectivamente. 


Ayuda en MATLAB 


Existen diversas alternativas para obtener ayuda sobre las numerosas funciones (y 
sus argumentos opcionales) implementadas en MATLAB. Para acceder a la ayuda de 
una función concreta de la que conocemos su sintaxis, podemos operar del modo si- 
guiente: o bien obtener la ayuda por pantalla (en la ventana de comandos), mediante 
la instrucción 


>> help nombre_funcion 
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o bien acceder a la misma a través de una ventana auxiliar, mediante la instrucción 
>> helpwin nombre_funcion 


Para desplegar el menú de ayuda genérico de MATLAB, se puede teclear directamente 
helpwin. De esta forma, se accede a una ventana de navegación que incluye un listado de 
tópicos, organizados por librerías, acerca de los cuales MATLAB dispone de tutoriales 
de ayuda. 


4. Variables y funciones elementales 


A diferencia de otros lenguajes de programación clásicos, MATLAB no necesita 
declarar variables. En principio, todas las variables son reales y basta con hacer uso de 
ellas para que queden automáticamente declaradas. Si escribimos 3 + 5 en el prompt de 
MATLAB, el resultado mostrado por pantalla será: 


>> 3+5 
ans = 


8 


Nótese que ans (apócope de answer) es una variable en la que MATLAB almacena el 
resultado de la última operación efectuada. 

Los comandos who y whos proporcionan información sobre las variables declaradas 
a lo largo de una sesión. En particular, who indica qué variables se están utilizando y 
whos incluye además el tamaño y tipo de cada variable. Por su parte, la instrucción 
clear permite borrar el contenido de todas las variables definidas durante dicha sesión. 
Para eliminar el contenido específico de la variable a, teclearemos clear a. MATLAB 
presenta además algunas variables especiales: 


o pi es el número T; 


o eps es el e de la máquina, es decir, el número positivo más pequeño tal que 
1+€H4H 1 en la aritmética del ordenador (eps = 2.2204e-016); 


o realmin y realmax son los números reales positivos más pequeño y más grande, 
respectivamente, que se pueden utilizar (realmin = 2.2251e-308 y realmax = 
1.7977e+308); 


o inf denota infinito; 


o NaN (not-a-number) es una magnitud no numérica que se obtiene, por ejemplo, 
como resultado de indeterminaciones matemáticas del tipo 0/0 ó oo — oo. 


En general, MATLAB almacena las constantes como matrices de dimensión 1 x 1 
y trabaja tanto con números reales como con números complejos. De hecho, i y j, 
por defecto, son variables que almacenan el valor de la unidad imaginaria y—1. Las 
funciones conj, abs y angle calculan el conjugado, el módulo y el argumento (en 
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radianes) de un número complejo dado, respectivamente. Por ejemplo, conj(3+4i) da 
como resultado 3-41. 

Los operadores aritméticos básicos son +, =-, *, / y 7, y hacen referencia a las opera- 
ciones de suma, resta, multiplicación, división y potencia, respectivamente. Adicional- 
mente, existe el operador X, que denota la inversa de la división (e.g., 211 da como 
resultado 0.5). La precedencia de operadores aritméticos es la habitual; no obstante, 
tal precedencia puede modificarse haciendo uso de paréntesis. 

Algunas de las funciones elementales implementadas en MATLAB son abs, sqrt, 
exp, log, log10, log2 y sign, así como las funciones trigonométricas sin, cos, tan, 
asin, acos y atan. Obsérvese que, en estas últimas, el ángulo viene expresado en 
radianes. 


5. Vectores y matrices 


Los vectores y matrices se definen en MATLAB introduciendo sus elementos entre 
corchetes. Los elementos de cada fila están separados por una coma (,) o un espacio 
en blanco, mientras que se utiliza un punto y coma (;) para separar dos filas consecu- 
tivas. Los siguientes comandos permiten definir un vector fila u de tres componentes, 
un vector columna v de cuatro componentes y una matriz A con tres filas y dos columnas: 


>u=[20-1]; 
> v = [3;-4;1;3]; 
[1 -1;4 3;2 10]; 


> A 


Alternativamente, los elementos de vectores y matrices pueden definirse a partir de uno 
o varios bucles for anidados, cuya sintaxis se especifica en la Sección 7. 

Por otra parte, MATLAB cuenta con una serie de comandos especialmente diseñados 
para la definición de vectores y matrices. Entre ellos, cabe destacar los siguientes: 


o a:incr:b genera un vector cuyos elementos son a, a+ incr, a+ 2x* incr, ..., 
a+kxincr, siendo k el mayor entero que satisface a+ k x* incr < b; en parti- 
cular, a:b considera incr = 1; 


o linspace(a,b,n) genera un vector de n componentes equiespaciadas, cuyo primer 
y último elementos son a y b, respectivamente; 


o rand(m,n) genera una matriz m x n, cuyos elementos son números aleatorios dis- 
tribuidos uniformemente entre O y 1; en particular, rand(n) genera una matriz 
cuadrada de orden n que contiene elementos de este tipo; 


o zeros(m,n) genera una matriz m x n, cuyos elementos son iguales a 0; en parti- 
cular, zeros(n) genera una matriz cuadrada de orden n formada por elementos 
nulos; 


o ones(m,n) genera una matriz m Xx n, cuyos elementos son iguales a 1; en particular, 
ones (n) genera una matriz cuadrada de orden n formada por elementos unidad; 
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o eye(m,n) genera una matriz m x n, cuyos elementos diagonales son iguales a 1 y 
cuyos elementos extradiagonales son iguales a 0; en particular, eye(n) genera la 
matriz identidad de orden n; 


o diag(v,k), con v vector de n componentes, genera una matriz cuadrada de orden 
n + |k|, cuya k-ésima diagonal contiene los elementos de v; el valor k = 0 hace 
referencia a la diagonal principal, mientras que los valores k > 0 y k < O denotan 
la k-ésima diagonal por encima o por debajo de la principal, respectivamente; en 
particular, diag(v) equivale a diag(v,0); 


o diag(A,k), con A matriz, genera un vector columna formado por los elementos de 
la k-ésima diagonal de A; en particular, diag(A) equivale a diag(A,0) y devuelve 
la diagonal principal de A. 


Una vez definida una matriz, el acceso a sus elementos se realiza mediante el uso de 
paréntesis. Por ejemplo, para hacer que el elemento ubicado en la fila ¿ = 2 y columna 
3 =2 de la matriz A sea igual a 8, escribiremos: 


>> A(2,2) = 8; 


Es importante tener en cuenta que el tamaño de una matriz se puede ir modificando 
dinámicamente, tal y como se muestra en la siguiente secuencia de comandos: 


>> A = ones(2,2); A(1,2) = -3; A(4,3) = -5 


A = 
1. -3 0 
1 1 0 
0 0 0 
0 0. -5 


Los dos puntos (:) se utilizan para acceder a partes determinadas de una matriz (filas, 
columnas o bloques). Por ejemplo, dada la matriz A anterior: 


o A(1,:) proporciona la primera fila de A; 


o A(:,2) proporciona la segunda columna de A; 


O 


A(1:3,2:3) proporciona la siguiente submatriz: 


[ACE,2) (1,35 4(2,2) A(2,3)5 (3,2) 1(8,39)] 


o A(1,2:3) = 0 iguala a O los elementos A(1,2) y A(1,3); 


A(:,3) = [] elimina la tercera columna de A; 


O 


A= [A;[1 21] añade a A una nueva fila; 


O 


o A=[A [1;2;3;4;5]] añade a A una nueva columna; 
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o A([1 31,:) =AC[3 1],:) intercambia las filas 1 y 3 de A sin necesidad de uti- 
lizar variables auxiliares. 


Ejercicio 1. Sea la siguiente matriz cuadrada: 
201 
A= 132 
140 


(a) Construye la matriz que se obtiene yuxtaponiendo la matriz identidad de orden 3 
a la derecha de la matriz A. A continuación, suma a la segunda fila de la matriz 
obtenida la primera fila multiplicada por 3. Por último, intercambia las columnas 
2 y 3 de la matriz resultante. 


(b) Construye dos nuevas matrices cuyas filas sean, respectivamente, las filas 1 y 3 
de A y las columnas 1 y 3 de A. 


Ejercicio 2. En una sola instrucción de MATLAB, sustituye todos los valores de la 
diagonal principal de una matriz cuadrada por cero. 


Ejercicio 3. En una sola instrucción de MATLAB, sustituye todos los valores de la 
diagonal principal de una matriz cuadrada por los elementos de un vector dado. 


5.1. Aritmética de vectores y matrices 


Los operadores aritméticos +, = y * permiten realizar sumas, restas y productos de 
vectores y/o matrices, siempre que éstos tengan las dimensiones adecuadas. En este 
caso, el operador / actúa de la siguiente forma: dadas las matrices A y B, el cociente 
A/B es equivalente al producto de A por la inversa de B. Además, mediante el operador 
” se pueden calcular potencias de una matriz cuadrada. Como novedad, es posible efec- 
tuar Operaciones componente a componente, añadiendo un punto delante del operador 
aritmético correspondiente. Por ejemplo, dadas dos matrices A y B de las mismas dimen- 
siones, A.* B tiene como resultado una matriz C cuyo elemento (+, 7) es el producto del 
elemento (+, 7) de A por el elemento (i, j) de B. Por su parte, A.72 tiene como resultado 
una matriz D cuyo elemento (+, j) es el cuadrado del elemento (i, j) de A. 

Por otra parte, el operador Y (usualmente denominado backslash) sirve para resolver 
sistemas de ecuaciones lineales. Dado un sistema genérico de la forma Ax = b, donde 
A es la matriz de coeficientes, x representa el vector de incógnitas y b denota el vector 
de términos independientes, MATLAB puede resolverlo mediante la instrucción: 


>> x= AWb 


Si el sistema es compatible determinado, la solución se obtiene mediante una factori- 
zación LU con pivotaje parcial (también denominada factorización PA = LU). En el 
caso particular de que A sea una matriz hermitiana y definida positiva, la factorización 
anterior es reemplazada por el método de Cholesky. Si el sistema es compatible inde- 
terminado, MATLAB devuelve una de sus infinitas soluciones. Por último, cuando el 


Polinomios el 


sistema es incompatible, se obtiene la solución que minimiza norm(A*x-b) aplicando 
el método de mínimos cuadrados?. 
Finalmente, dada la matriz A y los vectores u y v, se tiene: 


o [m,n] = size(A) devuelve las dimensiones de A; en particular, m indica el número 
de filas de A y n, su número de columnas; 


o rank(A) calcula el rango de A; 

o det(A) calcula el determinante de A, siempre que ésta sea cuadrada; 

o inv(A) calcula la matriz inversa de A, siempre que ésta sea regular; 

o trace(A) calcula la traza (suma de los elementos de la diagonal principal) de A; 


o [vec,val] = eig(A) devuelve una matriz llena (vec), cuyas columnas son los 
vectores propios de A, y una matriz diagonal (val), cuyos elementos diagonales 
son los valores propios de A, siempre que ésta sea cuadrada; 


o A” calcula la matriz traspuesta de A, si ésta es real, y su matriz traspuesta y 
conjugada, si es compleja; para calcular la matriz traspuesta sin conjugar de una 
matriz compleja A, se utiliza el comando A.?; 


o length(u) devuelve el número de elementos del vector u; si se aplica sobre la 
matriz A, devuelve el máximo entre su número de filas y su número de columnas; 


o dot(u,v) calcula el producto escalar de u y v; 


o cross(u,v) calcula el producto vectorial de u y v, siendo u y v vectores de tres 
componentes. 


Es interesante mencionar que las funciones elementales implementadas en MATLAB 
operan de forma natural sobre vectores y matrices: el resultado obtenido es el mismo 
que si se aplica la función elemento a elemento. Por ejemplo: 


>> a = [-pi/4 0; O pi/4]; taní(a) 


ans = 
-1.0000 0 
o 1.0000 


6. Polinomios 


En MATLAB, un polinomio se representa mediante un vector fila, cuyos elemen- 
tos son los coeficientes del polinomio en orden decreciente de grado (empezando por 
el coeficiente director y terminando por el término independiente). Así, el polinomio 
p(x) =2x* + 5x — 1 se representa mediante el vector fila p = [2 O 5 -11. Nótese que, 


3Dado un vector v, el comando norm(v) calcula la norma euclídea de dicho vector. 
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cuando algún coeficiente del polinomio es nulo, se debe asignar el valor O a la posición 
correspondiente del vector. 

Veamos a continuación algunas funciones de MATLAB específicas para polinomios. 
Dados un vector v, una constante k y dos polinomios p(x) y q(x), representados en 
MATLAB mediante los vectores p y q, respectivamente, se tiene: 


o polyval(p,k) evalúa p(x) en k; 


o polyval(p,v) devuelve, en un vector de la misma dimensión que v, las evalua- 
ciones de p(x) en cada una de las componentes de dicho vector; 


o roots(p) devuelve las raíces (reales y complejas) del polinomio p(1); 


o poly(vw) devuelve los coeficientes del polinomio mónico cuyas raíces son las com- 
ponentes de v; 


o polyder(p) devuelve los coeficientes del polinomio p'(1); 


o polyint(p,k) devuelve los coeficientes del polinomio f p(x) dx, con constante 
de integración (o término independiente) igual a k; en particular, polyint (p) 
considera una constante de integración igual a 0; 


o conv(p,q) devuelve los coeficientes del polinomio producto obtenido al multi- 
plicar p(x) por q(u); 


o [c,r] = deconv(p,q) devuelve los coeficientes de los polinomios cociente (vector 
c) y resto (vector r) obtenidos al dividir p(x) entre q(x1). 


Ejercicio 4. Dados los polinomios p(x) = 31*—1?+=x-—3 y q(x) = 4x* + 3a? — 
x +1, exprésalos según la sintaxis de MATLAB y súmalos utilizando las instrucciones 
necesarias. 


6.1. Ajuste por mínimos cuadrados 


El comando polyfit determina, mediante el método de mínimos cuadrados, el poli- 
nomio que aproxima a un conjunto de datos*. El grado del polinomio es indicado por 
el usuario. En concreto, si x e y son dos vectores de longitud m, p = polyfit(x,y,n) 
devuelve en p los coeficientes del polinomio de grado n que aproxima a los puntos dados 
por xe y. Sin =m- 1 y los elementos del vector x son distintos dos a dos, el polinomio 
obtenido pasa por todos los puntos? y recibe el nombre de polinomio interpolador de 
Lagrange. En el caso de que n >m-1,osin =m- 1 y el vector x contiene dos elementos 
iguales, el polinomio no es único y MATLAB devuelve un mensaje de aviso (warning). 

A modo de ejemplo, las siguientes sentencias calculan, mediante el método de míni- 
mos cuadrados, los polinomios de grado 2, 3 y 4 que aproximan a la colección de puntos 


%Es decir, dada una colección de puntos [(2;,y;)P”,, polyfit obtiene un polinomio p(x) de grado 
n tal que Y, (p(x;) — y;)? es mínimo. 
5Se dice que p(x) pasa por los puntos ((x;,y;))_, si verifica que p(1,) = y;, para i=1,2,...,m. 


Polinomios dd, 


FIGURA 1. Polinomios p2, p3 y p4 (de grado 2, 3 y 4, respectivamente) que aproximan a los 
datos (0,1), (3,4), (5,3), (7,5), (8,7)). Los datos se hallan representados por asteriscos. 


1(0,1),(8.4).(5, 3 (1,0118, 7) 
>x= [0357 81l; y = [1435 7]; 


>> p2 = polyfit(x,y,2) 
p2 = 
0.0322 0.3728 1.3394 
>> p3 = polyfit(x,y,3) 
pS = 
0.0696 -0.8010 2.7202 1.0212 
>> p4 = polyfit(x,y,4) 
p! = 


-0.0131 0.2929 -2.0012 4.7214 1.0000 


La Figura 1 muestra los datos del problema (representados por asteriscos) y los poli- 
nomios obtenidos. Obsérvese que el polinomio de grado 4 es el polinomio interpolador 
de Lagrange asociado a los cinco puntos dados. 

Debido al carácter fuertemente oscilante de los polinomios de orden alto, se suelen 
tomar valores bajos de n, aunque m sea grande. Por otra parte, si estamos interesados 
en realizar una interpolación polinómica que involucre a un número elevado de puntos, 
es más adecuado hacer uso del comando spline, que permite calcular splines cúbicos 
a partir de los datos del problema. Asimismo, existen comandos más generales para 
realizar interpolación en una, dos y tres dimensiones, denominados respectivamente 
interp1, interp2 e interp3, cuyas principales características pueden consultarse me- 
diante los comandos de ayuda de MATLAB. 
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7. Bucles y estructuras de decisión 


En MATLAB, la estructura: 
for indice = inicio:incr:final 
sentencias 


end 


implementa un bucle for, en el que las sentencias son ejecutadas repetidamente mien- 
tras la variable indice toma los valores del vector inicio: incr:final. Por su parte, 
la estructura: 


while condición 
sentencias 


end 


implementa un bucle while, en el que las sentencias son ejecutadas repetidamente 
mientras la condición es cierta. La estructura de decisión más importante en MATLAB 
es el comando if. Su sintaxis es la siguiente: 


if condición_1 
sentencias_1 
elseif condición 2 


sentencias_2 


else 
sentencias 


end 


Si la condición_1 es verdadera, se ejecutan las sentencias_1; en caso contrario, si la 
condición_2 es verdadera, se ejecutan las sentencias_2. Y se opera así sucesivamente. 
En caso de que ninguna de las condiciones indicadas se verifique, se ejecutan las sen- 
tencias asociadas al comando else. 


7.1. Operadores relacionales y lógicos 


Las condiciones que aparecen en los comandos while e if descritos en el epígrafe 
anterior suelen involucrar a los operadores relacionales y lógicos elementales. En MAT- 
LAB, tales operadores se representan del modo siguiente: 


mayor que > | mayor que o igual a >= | igual a == 
menor que < | menor que o igual a <= | distinto de «== 
y E 0 | [no o 
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En general, el resultado de una comparación es una variable lógica o booleana, cuyo 
valor es 1 si la comparación es verdadera y 0 si ésta es falsal, 


8. Ficheros script y funciones 


Hasta el momento, hemos descrito cómo trabajar en MATLAB de forma interactiva 
a través de la ventana de comandos. Sin embargo, también es posible (y muchas veces 
recomendable) utilizar los denominados M-files: ficheros de texto con extensión .m que 
contienen instrucciones en el lenguaje de programación propio de MATLAB. Existen 
dos tipos fundamentales de M-files: 


o ficheros script: contienen una secuencia de comandos de MATLAB que se ejecutan 
como si estuviésemos trabajando en modo interactivo; 


o funciones: son el equivalente en MATLAB a las subrutinas de los lenguajes de 
programación clásicos; toda función se declara mediante el comando function, 
que aparece siempre en su primera línea ejecutable. 


Para crear M-files, es posible teclear directamente edit nombre fichero en el prompt 
de MATLAB. Al ejecutar esta instrucción, se abre el editor de MATLAB y se crea el 
fichero nombre_fichero.m. Alternativamente, se pueden seleccionar las opciones New 
> Seript o New —> Function en la pestaña Home. En general, es conveniente guardar los 
M-files en la carpeta seleccionada como directorio actual de trabajo o Current Folder. 


8.1. Ficheros script 


Un fichero script es un M-file que agrupa una serie de instrucciones de MATLAB. 
La ejecución de un script no requiere argumentos de entrada ni proporciona argumen- 
tos de salida. Este tipo de ficheros está especialmente indicado para aquellos casos en 
los que necesitamos ejecutar un conjunto de órdenes repetidamente y queremos evitar 
teclearlas en cada ejecución. 


Ejercicio 5. Crea un fichero script que muestre por pantalla los valores de las variables 
eps, realmin y realmax”. 


8.2. Funciones 


MATLAB nos permite implementar nuestras propias funciones a través de ficheros 
con extensión .m de tipo función. En este caso, es recomendable asignar el mismo 
nombre al M-file y a la función en él implementada. 

La sintaxis de una función es la siguiente: 


SEl carácter - se puede obtener con las teclas A1t + 126 del teclado numérico. 
"Para ello, puedes utilizar el comando disp que encontrarás en la ayuda (véase también la Sección 
9 de este documento). 
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function [output1,output2,...] = nombre funcion(inputl1,input2,...) 


% La cabecera de la función contiene los comentarios que aparecen 
% al teclear help nombre_funcion en el prompt de MATLAB. 


Obsérvese que en la primera línea se especifica que el M-file en cuestión contiene a una 
función. En particular, se dice cuál es el nombre de la función y cuáles son sus argu- 
mentos de entrada (entre paréntesis) y de salida (entre corchetes). Las variables que se 
definen dentro del código de la función son locales y desaparecen una vez que termina 
su ejecución. Del mismo modo, la función no puede acceder a ninguna variable definida 
fuera de ella misma. El carácter % se utiliza para introducir comentarios. Finalmente, si 
una instrucción contiene un cambio de línea, deberemos escribir tres puntos suspensivos 
(...) antes de pasar a la línea siguiente. 


Ejercicio 6. Implementa una función denominada suma_vector, que admita un vector 
dado como argumento de entrada y devuelva la suma de sus componentes como argu- 
mento de salida. 


Ejercicio 7. El método de Gauss es el algoritmo más sencillo para la resolución de 
sistemas de ecuaciones lineales. En forma de pseudocódigo, puede expresarse del modo 
siguiente: 
fori=1:n-—1 
fork=:1i+1l:n 
L= Us) Qu 
forj=:1i+1:2n 
Akj —= QAkj — 2. Qs; 
end 
bj = bx —£ - dy 
end 
end 
Un = Da han 
fori=n-=1:-1:1 
n 
Ly = (a = y 042) /asi 


j=i+1 
end 


Implementa en MATLAB el método de Gauss mediante una función que presente la 
siguiente cabecera: 
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/ GAUSS Método de Gauss. 

% 

/ X = GAUSS(A,B) aplica el método de Gauss para resolver 
% el sistema A*X = B. 


/ Argumentos de entrada: 
% A: matriz de coeficientes; 
% B: término independiente. 


7 Argumento de salida: 
% X: solución. 


9. Entrada y salida de datos 


En esta sección, trataremos de mostrar las diversas alternativas que ofrece MATLAB 
para la adquisición de datos por parte del usuario y su salida por pantalla, así como 
para la lectura y escritura en ficheros de distinto tipo. 

El comando input permite imprimir un mensaje en la línea de comandos de MAT- 
LAB y obtener un valor numérico o el resultado de una expresión introducida por el 
usuario. Al ejecutar la instrucción, MATLAB imprime el mensaje por pantalla y espera 
a que el usuario teclee el valor numérico o la expresión. En el caso de que se introduzca 
una expresión, ésta es evaluada con los valores actuales de las variables de MATLAB 
(contenidas en el espacio de trabajo o workspace). Veamos un ejemplo de uso de esta 
función: 


>> n = input("Teclee el número de ecuaciones: ?) 
Teclee el número de ecuaciones: 10 
n= 

10 


Si a una instrucción de este tipo se le añade el argumento ”s?, el texto tecleado se lee 
y almacena en un string (o cadena de caracteres) sin evaluar expresiones. 

El comando disp permite mostrar una cadena de caracteres o el valor de una varia- 
ble en el formato actual. Para escribir en una misma línea texto y variables numéricas, 
se puede hacer uso del comando num2str, que convierte el valor de la variable especifi- 
cada en un string. Por ejemplo: 


>> disp([*El sistema tiene ” num2str(n) ” ecuaciones?]) 
El sistema tiene 10 ecuaciones 


Por otra parte, MATLAB cuenta con un asistente de importación de datos que 
nos permite importar datos de una forma sencilla, sin tener que realizar suposiciones 
previas sobre el formato de los mismos. Este asistente se activa seleccionando la opción 
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Import Data en la pestaña Home o tecleando la instrucción uiimport en la ventana de 
comandos. 

Veamos a continuación algunas funciones diseñadas para la lectura y escritura en 
ficheros ASCIT: 


o [fid,texto] = fopen(*nombre fichero?,”c”) abre el fichero nombre fichero, 
en el que se podrán llevar a cabo las operaciones determinadas por c, devuelve un 
identificador de fichero fid y, en el caso de producirse algún error en el proceso, 
un mensaje de error adicional en la variable texto; c puede tomar los valores 
r (lectura), w (escritura que reemplaza los contenidos previos), a (escritura que 
se añade a los contenidos previos) y sus correspondientes r+, w+ y a+ (lectura y 
escritura); 


o st = fclose(fid) cierra el fichero indicado por el identificador fid; st representa 
un valor de retorno para posibles errores; 


o [var1,var2,...] = fscanf(fid,” formato”) lee las variables var1, var2, ..., 
del fichero indicado por el identificador fid; la cadena formato contiene los es- 
pecificadores de formato para las variables, que pueden ser %s (cadena de carac- 
teres), %d (variable entera) o %f (variable de coma flotante), entre otros; 


o fprintf(fid, formato”, var1,var2,...) almacena la salida formateada en el 
fichero indicado por el identificador fid; la cadena formato contiene los especi- 
ficadores de formato de escritura, que son similares a los del lenguaje C; si el 
argumento fid se omite o tiene asignado un valor de 1, el comando fprintf 
puede utilizarse para escribir en la ventana de comandos. 


9.1. Intercambio de datos con Microsoft Excel 


Finalmente, es posible importar y exportar datos de otras aplicaciones, tales como 
Microsoft Excel98. En este contexto, instrucciones del tipo: 


>> x1sread(*nombre_fichero”,”nombre_hoja?,*rango?) 


permiten a MATLAB leer aquellos datos ubicados en la región determinada por rango 
de la hoja nombre_hoja incluida en el fichero de Excel nombre_fichero. 

Análogamente, para exportar datos de MATLAB a Excel, contamos con instruc- 
ciones del tipo: 


>> x1swrite(*nombre fichero”, nombre variable,'*nombre hoja”,*rango?) 


que permiten a MATLAB volcar el contenido de la variable nombre variable (que, en 
general, será una matriz) en la región determinada por rango de la hoja nombre_hoja 
incluida en el fichero de Excel nombre_fichero. 


SMicrosoft Excel es una marca registrada de Microsoft Corporation. 
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10. Symbolic Math Toolbox 


El paquete Symbolic Math Toolbox es uno de los muchos paquetes que amplían 
las funcionalidades de MATLAB”. Este paquete está basado en el núcleo del software 
Maple91%, que permite realizar todos los cálculos de forma simbólica y en precisión 
variable. Para acceder a las funciones contenidas en este paquete, basta con teclear 
help symbolic. 

En este contexto, los comandos sym y syms sirven para declarar objetos simbólicos. 
Así, con las siguientes instrucciones: 


>> syms X y; syms z real; sqroot2 = sym(*sqrt(2)?”); 


declaramos, respectivamente, dos variables simbólicas x e y, una variable simbólica z 

real y una constante simbólica sqroot2, cuyo valor es y2. Si se desea evaluar una 

constante simbólica const, es posible hacer uso de las instrucciones double(const) y 

vpa(const,d), donde d indica el número de dígitos de precisión requerido. Por ejemplo: 
>> format long; double(sqroot2) 


ans = 


1.414213562373095 
>> vpa(sqroot2,40) 
ans = 
1.41421356237309504880168872420969807857 


En cambio, para evaluar una expresión expr que contenga variables simbólicas (e.8., x 
e y), se utiliza la instrucción: 


>> subs (expr, 1x,y),1x0,y0)) 


Este comando evalúa la expresión expr mediante la asignación de los valores x0 e yO 

a las variables x e y, respectivamente. Existen instrucciones adicionales para llevar a 

cabo manipulaciones de expresiones simbólicas, tales como expand, simplify, factor 

y collect. Su definición y aplicaciones puede consultarse en la ayuda de MATLAB. 
Por su parte, la instrucción: 


>> solve(*eqn1”,*eqn2”,...,”eqnN”,*var1*,*'var2*,...,*varN?) 
permite resolver una ecuación o un sistema de ecuaciones, siendo eqn1, eqn2, ..., eqnN 
las ecuaciones que se desean resolver y var1, var2,..., varWN las incógnitas correspon- 


dientes. Veamos los siguientes ejemplos: 


>> syms x; solve(*exp(2*x)+3xexp(x)=54?) 


“Para conocer los paquetes contenidos en una instalación de MATLAB, basta con teclear ver en la 
ventana de comandos. 
l0Maple es una marca registrada de Waterloo Maple, Inc. 
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ans = 


log (6) 
pixi+log(9) 


>> syms Xx y; [x,y] = solve(*x"2+x*y+y=a*,*x" 2-4*x+3=0*,*x”,*y?) 


SS 
1 
3 

y = 
a/2-1/2 
a/4-9/4 


En la Sección 6, hemos visto cómo trabajar con polinomios a través de vectores. 
No obstante, MATLAB también puede operar con polinomios de forma simbólica. Los 
siguientes comandos convierten la expresión de un polinomio en formato vector a su 
expresión simbólica y viceversa: 


o poly2sym(u) devuelve la expresión simbólica del polinomio cuyos coeficientes 
vienen determinados por los elementos del vector fila u; 


o sym2poly(expr) devuelve un vector fila que contiene los coeficientes del polinomio 
cuya expresión simbólica viene dada por expr. 


Veamos a continuación una serie de instrucciones que nos van a permitir realizar diversas 
operaciones sobre expresiones simbólicas!!: 


o diff (expr,var,n) calcula la derivada n-ésima de la expresión expr con respecto 
a la variable var; el último argumento n es opcional, siendo 1 su valor por defecto; 


o int(expr,var) calcula la integral indefinida de la expresión expr con respecto 
a la variable var; análogamente, int(expr,var,a,b) calcula la integral definida 
de la expresión expr con respecto a la variable var en el intervalo [a, b); 


o limit(expr,var,a) calcula el límite de la expresión expr cuando la variable 
var tiende a a, siendo O el valor por defecto de a; para obtener límites laterales, 
teclearemos limit(expr,var,a,*dir”), donde dir será right o left, respecti- 
vamente, si el límite lateral se calcula por la derecha o por la izquierda; 


o symsum(expr,var,a,b) calcula la suma de la serie cuyo término general viene 
dado por expr, para los valores de la variable var entre a y b; 


o taylor(f,n,a) calcula el polinomio de Taylor de grado n — 1 de la función f, 
evaluado en el punto a; el valor por defecto de a es 0; 


11 En tales instrucciones, el argumento var denota la variable independiente considerada y es siempre 
opcional, siendo x su valor por defecto. 
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o fourier(f) calcula la transformada de Fourier*? de una función f que, por de- 
fecto, tiene como variable independiente x; el resultado obtenido es, por defecto, 
una función que depende de w. El cálculo de la transformada inversa de Fourier 
puede llevarse a cabo mediante el comando ifourier; 


o laplace(f) calcula la transformada de Laplace*? de una función f que, por de- 
fecto, tiene como variable independiente t; el resultado obtenido es, por defecto, 
una función que depende de s. El cálculo de la transformada inversa de Laplace 
puede llevarse a cabo mediante el comando ilaplace. 


Finalmente, mediante instrucciones del tipo dsolve(*eqn1”,*eqn2”,...,”eqnN?), 
se pueden resolver ecuaciones diferenciales o sistemas de ecuaciones diferenciales. En 
este caso, utilizaremos la letra mayúscula D para referirnos a la primera derivada, la 
combinación D2 para referirnos a la segunda derivada, y así sucesivamente. Las cons- 
tantes arbitrarias de la solución, cuando aparezcan, se denotarán por C1, C2, etc. Veamos 
algunos ejemplos: 


>> dsolve(*Dy-3xy=12”) 
ans = 


C1xrexp(3*t)-4 


>> [x,y] = dsolve(*Dx=x+4x*y” ,*Dy=-4*x+y?) 
Xx = 


C3xcos (4*t)*exp(t)+C2*sin(4*t)*exp(t) 


y = 
C2x*cos (4*t)*exp(t)-C3*sin(4*t)*exp(t) 


Obsérvese que no es necesario declarar la variable y como simbólica antes de utilizar 
dsolve. Asimismo, también es posible especificar las condiciones iniciales o de contorno 
que permitan fijar las constantes de la solución, utilizando la sintaxis 
dsolve(*eqn”,*cond1”,*”cond2”,...). 

Por último, la instrucción ezplot(expr, [xmin,xmax]) nos permite representar 
expr en el intervalo cerrado de extremos xmin y xmax; por defecto, dicho intervalo 
es |-27, 271]. 


Ejercicio 8. La siguiente ecuación describe la evolución de la temperatura T(t) de 
cierto objeto sumergido en un baño líquido a temperatura Ti(t): 


aT(t) - 
10 7 + T(t) NN To(t) 


Supongamos que la temperatura del objeto es, inicialmente, T(0) = 70%F, mientras 


12Dada una función f(x), su transformada de Fourier viene dada por F(w) = [%_ f(x) e ““?dz. 


3Dada una función f(t), su transformada de Laplace viene dada por L(s) = f," f(t) e=*dt. 
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que la temperatura del baño es de 170F. Utiliza MATLAB para resolver las siguientes 
cuestiones: 


(a) Determina la temperatura del objeto T(t). 
(b) ¿Cuánto tiempo transcurrirá hasta que T(t) alcance un valor de 168 *F ? 


(c) Representa T'(t) como una función del tiempo. 


11. Salidas gráficas 


En esta sección, se describen las herramientas disponibles en MATLAB para realizar 
representaciones gráficas en dos y tres dimensiones. Teniendo en cuenta la estructura del 
programa, las salidas gráficas se pueden interpretar como una representación adecuada 
de vectores y matrices (considerando a estas últimas como conjuntos de vectores fila o 
columna). MATLAB utiliza un tipo especial de ventanas para realizar sus operaciones 
gráficas; tales ventanas se denominan genéricamente ventanas gráficas y aparecen repre- 
sentadas por la etiqueta Figure N, donde N denota el índice asociado a una determinada 
ventana. 

La generación y manipulación de gráficas en MATLAB suele realizarse desde la 
ventana de comandos o a partir de ficheros script. Alternativamente, es posible utilizar 
de forma interactiva el editor gráfico; para más información sobre su uso, véase help 
plotedit o el menú correspondiente de la ventana gráfica (Tools => Edit Plot). 


11.1. Representaciones bidimensionales 


El comando plot sirve para representar gráficas de funciones y curvas en el plano. 
Dados dos vectores x e y de dimensión n, la instrucción plot (x,y) inicializa una ven- 
tana gráfica y representa en el plano las componentes del vector y frente a las de x. 
La idea que subyace a este proceso es la siguiente: para i=1,2,...,n, MATLAB une 
secuencialmente los puntos del plano (x(i),y(1)) mediante segmentos de recta, de 
forma que, para un número suficientemente elevado de puntos, la gráfica obtenida tiene 
apariencia suave. Por ejemplo, si queremos representar la función f(x) = xsen(x) en el 
intervalo [0, 917], teclearemos las siguientes instrucciones: 


>> x = linspace(0,9*pi,200); y = x.*sin(x); plot(x,y) 


obteniendo la gráfica que se muestra en la Figura 2. Obsérvese que la operación producto 
se realiza componente a componente sobre el vector x (i.e., el símbolo * está precedido 
por un punto). Otra opción para obtener una gráfica similar es hacer uso de la instruc- 
ción: 

>> fplot("x*sin(x)”,[0,9*pi]) 


Finalmente, también es posible declarar: 


>> ezplot("x*sin(x)”, [0,9*pi]) 
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FIGURA 2. Representación gráfica de la función f(x) = xsen(z), en el intervalo [0, 97]. 


El comando plot permite asimismo generar gráficas de curvas definidas paramétri- 
camente, €.8.: 


>> t = 0:.001:2*pi; x = cos(3*t); y = sin(2*t); plot(x,y) 


genera la gráfica mostrada en la Figura 3. Obsérvese que, en este caso, aunque los valores 
de las variables x e y pertenezcan al mismo intervalo [—1, 1], las marcas mostradas en 
los ejes de abscisas y ordenadas son diferentes. En la siguiente subsección veremos cómo 
modificar esta propiedad de los ejes de coordenadas. 

Las representaciones gráficas obtenidas con el comando plot pueden personalizarse. 
La manera más sencilla de hacerlo consiste en añadir a dicho comando un tercer argu- 
mento, de modo que éste adquiere la forma plot (x,y,”linea”). En este caso, linea 
es una cadena de caracteres que especifica el formato de línea, el tipo de marca que se 
coloca sobre los puntos (x(1),y(1)) y/o el color utilizado. A continuación se detalla 
la sintaxis de estos elementos: 


o formato de línea: sólido (-), discontinuo (--), punteado (:), raya-punto (-.); 


o tipo de marca: punto (.), cruz (+), asterisco (*), aspa (x), círculo (o), cuadrado 
(s), diamante (d); 


o color: amarillo (y), magenta (m), rojo (r), verde (g), azul (b), cyan (c), 
blanco (w), negro (k). 


Por defecto, cada llamada a la función plot genera una gráfica que elimina la an- 
terior, en el caso de que exista. Para superponer varias gráficas en una misma figu- 
ra, se puede hacer uso de la instrucción hold on. El comando hold off desactiva 
tal funcionalidad. La Figura 4 muestra simultáneamente las gráficas de las funciones 


—p2 
fa) =e* y g(a) = e > /2, que han sido obtenidas utilizando el comando hold on. 


Más adelante incluiremos la secuencia de instrucciones completa que da lugar a dicha 
figura. 
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FIGURA 3. Representación gráfica de la curva (cos(3t),sen(2t)), para t € [0, 271]. 


Para manipular propiedades más específicas de las representaciones gráficas, es posi- 
ble incluir argumentos adicionales en el comando plot. La forma general de declarar 
dichos argumentos es la siguiente: 


>> plot(x,y,*linea”,”nombre1”,valor1,'*nombre2*,valor2,...) 


Algunos ejemplos de propiedades de este tipo son LineWidth, que controla el grosor de la 
línea, o MarkerSize, MarkerFaceColor y MarkerEdgeColor, que permiten especificar, 
respectivamente, el tamaño de las marcas, el color de su interior y el color de su borde. 
Una forma alternativa de indicar el color de la línea es añadir la propiedad Color, 
seguida de uno de los caracteres *y?, "m”, *r?, *g”, *b”, *c”, "w* o *k” especificados 
anteriormente, o de un vector de tres componentes (con valores entre O y 1) que define 
el color según el estándar RGB. 

Finalmente, existe una serie de comandos más generales que permiten controlar el 
entorno gráfico. Entre los más elementales, cabe destacar los siguientes: 


o title(*texto?) define el título de la gráfica; 


o xlabel(*texto*) e ylabel(*texto”) asignan las etiquetas correspondientes a 
los ejes x e y, respectivamente; 


o axis([xmin xmax ymin ymax]) reescala el tamaño de los ejes x e y; 


o text(x,y,'texto”) sitúa el texto especificado en las coordenadas (x, y) de la 
ventana gráfica; 


o gtext(*texto”) sitúa el texto especificado en las coordenadas de la ventana 
eráfica que se indican utilizando el ratón; 


o axis on muestra los ejes x e y; axis off oculta ambos ejes; 


o grid on muestra una malla superpuesta a la gráfica; grid off oculta la malla 
anterior; 
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FIGURA 4. Representación gráfica de las funciones f(x) =e"" y g(1) =e 2 /2, 


en el intervalo [-2, 2]. 


o legend(*texto1”,*texto2*,...) define una leyenda para la gráfica; 
o whitebg(*color”) asigna un color de fondo a la representación. 


Volviendo a la Figura 4, la siguiente secuencia de instrucciones genera la representación 
eráfica contenida en la misma: 

x=-2:.2:2; 

plot (x,exp(-x."2),*--*k?) 

hold on 

plot (x, expl-2%."2/2)/2,*=0b”) 

axis([-2,2,0,1.2]) 

title(*”Dos exponenciales?) 

xlabel (*x”) 

ylabel (y?) 

grid on 

legend(*e”]-x2"2)?,%e "1-2 "2/2)/2*) 

whitebg([0.9 0.9 0.9]) 
Nótese que la ausencia del prompt de MATLAB en la lista de comandos previa nos 
indica que ésta se encuentra definida en un fichero script. 

La visualización de diferentes subventanas gráficas en una misma ventana puede lle- 
varse a cabo mediante el comando subplot. En particular, la instrucción subplot (425) 
o, equivalentemente, subplot (4,2,5), divide la ventana gráfica en una matriz de sub- 


ventanas de dimensiones 4 x 2 y accede a la quinta subventana, es decir, a aquélla 
situada en la posición (3,1) de la matriz. Si, a continuación, tecleamos subplot (428), 
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estaremos especificando el acceso a la octava subventana, ubicada en la posición (4, 2) 
de la matriz. En la Figura 7, se muestra un ejemplo de uso del comando subplot; el ac- 
ceso a cada una de las subventanas se realiza mediante las instrucciones subplot (211) 
y subplot (212). 

Para concluir esta subsección, cabe mencionar algunos comandos adicionales dispo- 
nibles en MATLAB para la elaboración de representaciones gráficas bidimensionales, 
tales como plotyy, polar, semilogx, semilogy, loglog, stem, stairs, bar, area, 
line, fill o patch, cuya sintaxis y aplicación pueden consultarse en la ayuda. 


Ejercicio 9. Representa en MATLAB la curva epicicloide definida paramétricamente 
por las ecuaciones: 


x(t) = (a + b) cos(t) — bcos (221), 
y(t) = (a +b)sen(t) —bsen(2%t), 


donde a = 12 yb = 5, para un valor del parámetro 0 <t < 107. Considera la misma 
escala para los ejes de abscisas y ordenadas e incorpora una malla adicional a la figura 
resultante. Por último, añade un título a la gráfica y especifica los nombres de ambos 
ejes. 


Ejercicio 10. Representa en una misma figura la gráfica de la función f(x) = e” y 
la de su polinomio de Taylor de segundo grado centrado en el origen. Identifica am- 
bas gráficas mediante la correspondiente leyenda y añade a la figura cuantos atributos 
estimes oportunos. 


11.2. Representaciones tridimensionales 


El comando plot3 es el análogo de plot en tres dimensiones. Si x, y y z son tres 
vectores de la misma dimensión n, la instrucción plot3(x,y,z) representa en el espacio 
tridimensional la curva que pasa por los puntos (x(i),y (4),z(1)), parai =1,2,...,n. 
Tales vectores suelen introducirse en forma paramétrica. Así, la secuencia de instruc- 
ciones: 


>>t = 0:0.1:10x*pi; 
x = exp(-t/20).*cos(t); y = exp(-t/20) .*sin(t); z = t; 
plot3(x;y,z) 


representa la espiral que se muestra en la Figura 5. El comando plot3 admite los 
mismos argumentos opcionales descritos para plot. Adicionalmente, mediante el botón 
Rotate 3D, ubicado en la barra de herramientas de la ventana gráfica, es posible realizar 
rotaciones de cualquier gráfica tridimensional. 

MATLAB es capaz de representar superficies en el espacio a través del coman- 
do surf. Para ilustrar su uso, vamos a representar la función definida como f(r) = 
sen(r)/r, donde r(x,y) = y 1? + y? +e, en el cuadrado [—10, 10] x [—10, 10]. Para ello, 


teclearemos: 


Salidas gráficas 25 
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FIGURA 5. Representación gráfica de la curva (0720 cos(t), e” Bsen(t), t), para t € [0, 107]. 


>> xx = -10:0.5:10; yy = xx; [lx,y] = meshgrid(xx,yy); 
r = sqrt(x."2+y."2)+eps; z = sin(r)./r; surf(x,y,z); 
colormap(*pink”); colorbar 


La salida gráfica obtenida aparece representada en la Figura 6. Nótese que e es un valor 
positivo próximo a O que evita que se anule el denominador de la función; en particular, 
utilizaremos la variable eps de MATLAB, cuyo valor es 2.2204e-016. En la secuencia 
de instrucciones previa, los vectores xx e yy definen sendas particiones sobre los lados 
del cuadrado. Por su parte, la función meshgrid genera dos matrices x e y: la matriz x 
presenta tantas filas como elementos contiene el vector yy, y almacena en cada una de 
ellas los elementos del vector xx; análogamente, la matriz y presenta tantas columnas 
como elementos contiene el vector xx, y almacena en cada una de ellas los elementos 
del vector yy. Para un caso sencillo, la idea es la siguiente: 


>> xx = [-1;0;11; yy = [-2;21; [x,y] = meshgrid(xx,yy) 


Xx = 
-1 0 
-1 0 1 

y = 
-2 -2 -2 
22 2 


El comando colorbar muestra una escala de color graduada, que se sitúa a la 
derecha de la figura e informa sobre el valor numérico correspondiente a cada color. 
Por su parte, el mapa de color utilizado puede modificarse mediante la instrucción 
colormap. MATLAB dispone de los siguientes mapas de color predefinidos: jet (por 
defecto), hot, cool, hsv, pink, copper, flag, gray y bone. Para un mapa de color 
específico, la distribución de color sobre una superficie se controla mediante el comando 
shading, que presenta las opciones shading faceted (por defecto), shading interp 
y shading flat. 

Como alternativa al comando surf, MATLAB dispone de la instrucción mesh para 
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FIGURA 6. Representación gráfica de la función f(r) =sen(r)/r, donde 
r(1,y) = yx? + y? +e, en el intervalo [-10, 10] x [-10, 10]; el 
parámetro e toma el valor de la variable eps de MATLAB. 


representar superficies en el espacio. Su forma de operar es muy similar a la descrita 
para surf y sirve para visualizar el mallado cartesiano que une los puntos que definen 
la superficie. 


Ejercicio 11. Representa en MATLAB la siguiente función: 


f(x, y) = sen(y? + x) — cos(y — a), 


en el intervalo 0O<x<T,0<y< rr. En particular, muestra en una misma ventana 
gráfica las salidas de los comandos surf y mesh, dispuestas horizontalmente una al la- 
do de la otra. Añade un título a cada una de las gráficas y especifica en cada caso los 
nombres de los ejes de coordenadas. 


La representación de las curvas de nivel de una función de dos variables se lleva a 
cabo en MATLAB mediante el comando contour. Una variante del mismo, dada por la 
instrucción contourf, permite colorear el espacio delimitado por las curvas represen- 
tadas. Para añadir texto sobre una determinada curva, es posible utilizar el comando 
clabel. Á continuación se muestra un ejemplo de aplicación de estas tres instrucciones. 
La salida gráfica correspondiente aparece representada en la Figura 7: 

xx = -2:0.1:2; 

yy = -2:0.1:2; 

[x,y] = meshgrid(xx, yy); 

subplot (211) 

c = contour(x,y,-3*y./(x."2+y."2+1),10); 

colormap(*pink?) 

colorbar 


clabel (c) 
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FIGURA 7. Curvas de nivel de la función f(x, y) = El en el intervalo 


[-2,2] x [-2,2]. Las gráficas superior e inferior corresponden a los 
comandos de MATLAB contour y contourf, respectivamente. 


subplot (212) 

ec = contourt(x,y, yx. 247. 2+1),10); 

colorbar 

clabel (c) 
El último argumento de las funciones contour y contourf permite fijar el número de 
curvas de nivel que se visualizan; se trata de un argumento opcional que, en este caso, 
toma un valor de 10. 

Finalmente, la representación de campos vectoriales en el plano y en el espacio se 


realiza a través de los comandos quiver y quiver3, respectivamente. En este contexto, 
la siguiente secuencia de instrucciones genera la salida gráfica de la Figura 8: 


xx = -1:.2:1; 

yy = -1:.2:1; 

[x,y] = meshgrid(xx, yy); 
quiver(x,y,-y/10,x/10,1.1) 

axis equal 

axis([-1 1-1 1) 

set(gca, *XTick*,-1:0.5:1) 

setígca, *ATickLabel*,?*-1*,*-0,5*,?0*,*0.5*,1>) 
set(gca, *YTick*,-1:0.5:1) 

selígea, Ylscklabel*,*=12,*-006%,70%,"0,6%,*17) 
title(*”Campo vectorial F(x,y)=(-y/10,x/10)?”) 
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Campo vectorial F(x,y) = (-y/10,x/10) 


e 


FIGURA 8. Representación gráfica del campo vectorial F(x,y) = (- 
en el intervalo [-1,1] x [-1, 1]. 


ale 
al* 
_— 


Cuando se representan campos vectoriales, MATLAB escala automáticamente los vec- 
tores con el fin de evitar que éstos se solapen entre sí. Esta propiedad puede modificarse 
mediante un último argumento opcional que, tanto en quiver como en quiver3, deter- 
mina el escalado deseado. Si dicho argumento es 0, el escalado automático se elimina, 
mostrándose la gráfica a escala 1: 1. 

Obsérvese que la instrucción set(gca,*XTick”,-1:0.5:1) sirve para localizar las 
marcas que aparecen a lo largo del eje de abscisas; en este caso, dichas marcas se sitúan 
en los extremos del intervalo [—1,1] y en los puntos interiores espaciados entre sí 0,5 
unidades. El argumento gca devuelve un identificador del último eje seleccionado. Por 
su parte, set(gca, *XTickLabel”,*-1*,*”-0.5*,*0”,*0.5*,*1”) añade las etiquetas 
correspondientes a los puntos mostrados. Los argumentos homólogos para el eje de 
ordenadas son YTick e YTickLabel. 

Por último, con el fin de ilustrar el uso del comando quiver3, vamos a representar 
la trayectoria de una partícula, cuya posición en función del tiempo (x(t), y(t), z(t)) 
viene dada por las expresiones x(t) = 2t, y(t) = 3t y z(t) = 10t — 16t?, para t € [0, 1]. 
La secuencia siguiente permite generar la gráfica representada en la Figura 9: 


t=.0:0.1:1; 
X = 2x*t; 
y = 3*t; 


z = 10*t-16*t .2; 
u = gradient(x); 
v = gradient(y); 
w = gradient(z); 
quiver3(x, y,Z,U,17,,0) 


view([120 40]) 
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FIGURA 9. Trayectoria de una partícula, cuya posición viene dada por 
(x(t), y(t), z(t)) = (2t, 3t, 10t — 16t?), para t e [0, 1]. 


En este caso, dado un vector v, el comando gradient (v) calcula la derivada numérica 
de las componentes del vector. 

Para concluir, cabe mencionar que MATLAB dispone de instrucciones más específi- 
cas para realizar representaciones tridimensionales. Entre ellas, destacan los comandos 
trimesh y trisurf, análogos a mesh y surf, que permiten representar superficies sobre 
dominios discretizados por mallados triangulares. 
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